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Nonspinning searches for spinning BH-NS binaries in ground-based detector data: Amplitude and 
mismatch predictions in the constant precession cone approximation 
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Current searches for compact binary mergers by ground-based gravitational-wave detectors assume for sim- 
plicity the two bodies are not spinning. If the binary contains compact objects with significant spin, then this 
can reduce the sensitivity of these searches, particularly for black hole-neutron star binaries. In this paper we 
investigate the effect of neglecting precession on the sensitivity of searches for spinning binaries using non- 
spinning waveform models. We demonstrate that in the sensitive band of Advanced LIGO, the angle between 
the binary's orbital angular momentum and its total angular momentum is approximately constant. Under this 
constant precession cone approximation, we show that the gravitational-wave phasing is modulated in two ways: 
a secular increase of the gravitational-wave phase due to precession and an oscillation around this secular in- 
crease. We show that this secular evolution occurs in precisely three ways, corresponding to physically different 
apparent evolutions of the binary's precession about the line of sight. We estimate the best possible fitting factor 
between any non-precessing template model and a single precessing signal, in the limit of a constant precession 
cone. Our closed form estimate of the fitting-factor depends only the geometry of the in-band precession cone; 
it does not depend explicitly on binary parameters, detector response, or details of either signal model. The 
precessing black hole-neutron star waveforms least accurately matched by nonspinning waveforms correspond 
to viewing geometries where the precession cone sweeps the orbital plane repeatedly across the line of sight, in 
an unfavorable polarization alignment. 
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I. INTRODUCTION 

Ground based gravitational-wave detector networks (no- 
tably LIGO lHJ-O and Virgo HI |5l) have performed several 
searches for the inspiral and merger of binaries containing 
black holes (BH) and neutron stars (NS) 16- 151 . The co- 
alesence of NS-NS, BH-NS, and BH-BH binaries are the 
most promising sources of gravitational-waves for these net- 
works |16|. For the lowest-mass compact binaries M = 
nil + ?Ti2 ^ 15Mq, the response of the detector to a bi- 
nary merger with arbitrary masses and spins is well under- 
stood IIT7H23JI . To identify the presence and properties of a 
signal, the classic approach to data analysis on nonspinning 
binaries has been matched filtering, linearly projecting the 
data against each member of a template bank: a discrete ar- 
ray of signal models, spaced to cover all possibilities with 
minimal inter-template and template-signal mismatch ll24l - 
l28l . For nonspinning binaries, the templates needed lie in 
the two-dimensional signal manifold itself (i.e., in rrii, m2 or 
more naturally in the chirp mass Mc — Mtfl^ and symmet- 
ric mass ratio 77 = toiTO2/M^) 1241 |25l l27l . For spinning bi- 
naries with more physical degrees of freedom, the added phys- 
ical parameters require either (i) exact template banks that are 
impractically large, with a significant increase in computa- 
tional burden and false alarm rate; or (ii) the use of an approx- 
imate detection template family, whose waveforms usually 
at best imperfectly approximate the target waveforms |18- 
l2n l29i but whose template banks are smaller and better un- 
derstood. Though the more generic BCV detection template 
family ll20l has been applied in previous searches for spin- 
ning systemsQO, 30 1, recent studies suggest that it is no more 
effective than the nonspinning template bank for low-mass 
single-spin binaries |26|. Current templated searches for the 



gravitational-wave signature of low-mass merging binaries are 
peiformed using this simple nonspinning bank fT2l [T3l [TSl . 
However, these results show that searching for generic spin- 
ning binaries with a nonspinning search is sub-optimal and 
further work is needed to construct an optimal search for spin- 
ning binaries in Advanced LIGO and Advanced Virgo. Such 
a search may involve, e.g. a hierarchical combination of a 
nonspinning search and a search using spinning templates, 
or a nonspinning search in some regions of the parameter 
space, with targeted spinning searches in others. To devise 
the correct stratergy, it is essential to understand where and 
why advanced-detector searches will lose signal-to-noise ra- 
tio when searching for spinning binaries using nonspinning 
waveforms; this is the goal of this paper 

Recent work on binaries where the spin is aligned with the 
orbital angular momentum has demonstrated that advanced 
detectors are more sensitive to the effects of spin on the or- 
bital phase than first generation detectors 131]. However, it 
has been shown that a slightly more generic search would re- 
cover all nonprecessing signals well by capturing the phase 
evolution of generic nonprecessing sources [31]. Nonethe- 
less, a noticeable fraction of these spinning binary mass, spin, 
and orientation parameters A are at best poorly matched with 
nonspinning and even generic nonprecessing templates, some 
with match lower than ~ 0.6. These large mismatches occur 
because nonspinning templates cannot capture the effects of 
spin. 

This paper is concerned with binaries where the spin and 
angular momentum of the binary are not aligned. This induces 
precession, causing the orbital angular momentum direction L 
changes with time. Because the instantaneous beam-pattern 
of the gravitational-wave emission is oriented along L, the 
changing angle between L and the line of sight n produces 



large amplitude and phase modulations that are absent in non- 
precessing signals. Poor matches with nonspinning templates 
are disproportionately concentrated in systems that exhibit 
precession: systems with significant spin-orbit misalignment 
and asymmetric mass ratio. Conversely, binaries with tight 
spin-orbit alignment or simply limited variation in L • n are 
well-matched by non-precessing templates that capture the ef- 
fects of aligned spin ' on the gravitational-wave phase ll32l - 
[34l . Understanding the effect of precession is key to de- 
veloping an optimal search for spinning binaries in second- 
generation gravitational- wave detectors. 

In this paper, we investigate the gravitational-wave sig- 
nal from precessing black hole-neutron star binaries. Except 
for a short epoch of transitional precession (typically when 
\L\ ~ \S\; see Il35j ). these binaries usually undergo simple 
precession, where the orbital angular momentum precesses 
around the total angular momentum with an ever-increasing 
opening angle 1 35 1 . By investigating the spin evolution equa- 
tions of Apostolatos et al. (351, we find that opening angle 
of the precession cone (i.e., the angle between the total angu- 
lar momentum J and the orbital angular momentum L) does 
not expand significantly while the binary's signal is in the de- 
tector's sensitive band. We perform analytic studies under the 
assumption that this opening angle is exactly constant; we call 
this the constant precession cone approximation. We find the 
gravitational-wave phasing is modulated in two ways. The 
first way that precession modifies the signal is via a secular 
change of the gravitational-wave phase. In addition to the ef- 
fect of orbital dynamics, the gravitational-wave signal also ac- 
cumulates a geometrical phase, a fraction of the "precession 
phase" of L about J. Previous work has demonstrated that 
precession induces a secular change in the phase accumulat- 
ing over each precession cycle, both in post-Newtonian stud- 
ies |fT9ll35l and in numerical simulations of GR |36|. Though 
long recognized, this factor is not included in standard an- 
alytic non-precessing template models If32|j34l .^ We further- 
more find this secular evolution occurs in three different ways, 
corresponding to physically different apparent evolutions of 
the precession about the line of sight. These three ways corre- 
spond precisely to whether the cone swept out by L encloses 
zero, one, or two of the "null lines", which are the directions 
that produce exactly zero instantaneous response in the detec- 
tor The three regions correspond to whether zero, one, or both 
directions are enclosed by the path swept out by i - hence- 
forth, the "precession cone." Each type of configuration leads 
to a unique effect on the phasing of the binary. If you are 
within one of the regions, you accumulate secular phase at the 
same rate. 



' Aligned spin means that tlie spins of the compact objects are aligned with 
the orbital angular momentum. In this case there is no precession. 

^ The full effects of precession are of course captured in full time-domain 
calculations of orbit dynamics and precession. These direct simulations, 
however, are presently too slow to be used in a search. Additionally, no 
one has constructed a complete template bank (a discrete set of reference 
signals) that adequately but minimally covers all possible signal options. 



The second way that precession modifies the signal is by in- 
troducing phase (and amplitude) modulations, on top of sec- 
ular evolution. These modulations cannot be captured by a 
non-precessing template. A non-precessing signal model will 
therefore not match a precessing signal optimally. Moreover, 
the highest possible mismatch between a precessing signal 
and non-precessing model can be estimated geometrically. 
In other words, our calculation tells us how easily a non- 
precessing search can find each particular precessing signal. 
For simplicity, we perform our analytic studies assuming a 
source directly overhead a single detector. Our results will be 
generalized to a full network in a future analysis. Though 
simplified, our study is nonetheless directly applicable to 
real searches.^ Moreover, given the complexity of realistic 
searches, special cases where searches have well-understood 
performance provide invaluable tools for code validation and 
search calibration. Past, present, and future searches therefore 
benefit from the goal of this paper: a simple model for how 
well nonspinning searches recover most spinning, precessing 
low-mass binaries. 

We will perform a detailed point-by-point numerical com- 
parison of this expression to a large-scale Monte Carlo study 
of synthetic searches in a subsequent publication. To compare 
with existing results in the literature, however, we apply our 
purely analytic results to randomly selected binary parame- 
ters. We find our expressions quantitatively agree with previ- 
ous Monte Carlo studies of nonprecessing searches 1.20, .37J , 
including recent investigations that include both mass and 
aligned spin as parameters fy\\. For instance, adopting the 
same spin distribution as 1311 , we reproduce a similar distri- 
bution of fitting factors. Unlike previous purely Monte Carlo 
studies, however, our analytic approach allows us to quanti- 
tatively predict precisely which systems are not well fit with 
nonprecessing templates. The performance of a nonprecess- 
ing search on BH-NS systems with any mass and spin distribu- 
tion can be easily predicted using our amplitude and mismatch 
predictions, and the simple predictions for the parameters that 
cause bad matches point the way toward improved searches. 

In Section [n] we demonstrate that BH-NS binaries seen by 
gravitational-wave detectors will often evolve on a constant 



precession cone. In Section III we parameterize the ways in 
which the binary's precession cone can be oriented relative to 
a single detector and its line of sight to the source. We find 
that there are three distinct regions in this configuration space. 



In Section IV we review the leading-order gravitational-wave 



emission from a precessing binary, decomposing it into sec- 
ular and modulated factors. We show how the gravitational- 
wave signal from each region has a distinctive secular depen- 
dence. Finally, in Section IV] we average the amplitude over 
the precession cone, relating the mean power seen along the 
line of sight to the power expected from an optimally ori- 
ented source. Additionally, we argue that non-precessing sig- 



' For example, the two-detector LIGO network is nearly aligned and approx- 
imately sensitive to one polarization at a time. In a subsequent paper we 
will address explicitly how to rescale our results to the two-detector LIGO 
network. 



nal templates necessarily cannot reproduce oscillations in the 
precessing signal model. Using a similar average and assum- 
ing a generic signal model reproduces the secular phase, we 
find an expression for the mismatch between non-precessing 
and precessing signals. 



n. ORBIT DYNAMICS AND THE CONSTANT 
PRECESSION CONE LIMIT 

We consider a binary of two compact objects with masses 
rrii and m2, with nii > TO2, and intrinsic spins 5*1 and 6*2. 
From the coordinate separation r and the velocity v ~ dtf, 
we can define the coordinate Newtonian angular momentum 
L = ufx V, where a ~ mim^ j^ jjjg reduced mass. 

^ ^ mi+m2 

For generic binaries, the orbit and spins evolve according to 
complicated position-, spin-, and velocity-dependent expres- 
sions Il32ll38l[39l . Averaging these post-Newtonian spin evo- 
lution equations over a few orbital periods leads to adiabatic 
expressions for the evolution of spin and orbital frequency 
[see, e.g., 135] (henceforth ACST), Schnittman ||39l 1: 
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In this expression we have omitted the quadrupole- 
monopole term described in ll39l . to more clearly correspond 
with the expressions from ACST on which we rely. We have 
included the spin-spin terms (containing both 5i and 5*2) for 
reference but we will soon assume that only one body is spin- 
ning. These evolution equations cause Si and ^2 to precess 
but leave their magnitude constant. The direction L of the 
orbital angular momentum also precesses and gravitational- 
wave emission causes the length of L to shrink as the binary 
evolves. 

As described in ACST, when either one mass is nonspin- 
ning (e.g., 6*2 — 0), the component masses are comparable 
(toi ~ 'rn2), or the smaller mass has a miniscule mass and 
necessarily spin {mi ^ 7712), these expressions reduce to ef- 
fectively single spin evolution, meaning equations of motion 
are equivalent to the dynamics of L and a single (effective) 
spin. For single spin binaries, spin-orbit evolution can occur 
in effectively two phases: (i) rare transitional precession when 
J = Z + S* ~ 0, the relative change in J is large, and J can 



change direction; and (ii) simple precession, when are nearly 
constant L and S precess roughly steadily around J. 

Single-spin evolution is particularly appropriate for BH-NS 
binaries. Neutron stars likely do not spin rapidly {S2 ~ 0) - 
observations suggest spin periods of no less than a few mil- 
liseconds |40|. Observations and theoretical considerations 
constrain plausible NS masses to 1 — SAf©, ensuring the sec- 
ond object's minimal spin has little impact on the orbit about 
an a priori higher-mass companion. Albeit with less accuracy, 
BH-BH systems also often approximately satisfy the single- 
spin condition, as binaries that are either comparable mass 
or dominated by a single spinning body occupy most of the 
TTii , 7712 , S'l , 5*2 parameter space; see, e.g., the discussion in 

m. 



A. Simple precession 

For the circular orbits expected in astrophysical scenarios, 
the spin evolution equations usually imply simple precession: 
the orbital and spin angular momenta precess about the total 
angular momentum, which is fixed, except for negligible pre- 
cession of its own about its average direction [ACST]. Assum- 
ing ^2 = 0, the system evolves with k = L ■ Si and 1 5*1 1 both 
constant. At early times, the total angular momentum J is 
dominated by and closely aligned with L; as radiation carries 
away orbital angular momentum, eventually J is dominated 
by and closely aligned with S. As the direction of the total 
angular momentum is nearly conserved, this process involves 
substantial change in the spin and orbital angular momentum 
directions. Following ACST, we will employ the ratio of Si to 
the Newtonian angular momentum L = jirx w as a parameter: 



7(0 ^ \Si\/\L{t)\ = 



Xi-ml 



r]My/Mr{t) 



(2) 



Using this parameter, the opening angle angle /3 of the preces- 
sion cone (denoted Xl in ACST) can be expressed trigono- 
metrically as 



m 



arccos J ■ L = arccos 



1 + K7 



^1 + 2k7 + 72 



(3) 



Note that for a spin-dominated binary (-f ^ 1) k — cos /?; for 
an angular momentum dominated binary (7 ~ 0) /3 ~ 0; and 
for a spin-orbit aligned binary (k = 1) all three of the spin, 
orbit, and total angular momenta are aligned (/3 = 0). At the 
level of accuracy of interest here,^ the opening angle /3(/) at a 
given gravitational frequency / can be estimated by the above 



■* As described at greater length in a subsequent publication, our constant pre- 
cession cone approximation is designed to be accurate to of order several 
percent in power and mismatch. Quantitative comparisons of our approxi- 
mation to numerical simulations will be provided in that publication. 




FIG. 1: Constant precession cone approximation works: Opening 

angle /3(/) versus frequency for a lOi\f0 + 1.4A/0 binary witli xi = 
1 and arccos k = 0, vr/lO, . . . (bottom curves) up to arccos k ~ tt 
(top, rapidly changing curve). We estimate /3(/) using Eqs. j3|4| . 
This plot shows that except for highly misaligned binaries (arccos k 
large), /?(/) is nearly constant, changing at most a fraction of a radian 
across the sensitive band of present and future detectors. In general, 
only a small subset of masses and highly misaligned spins are finely 
tuned enough to produce significant /3 evolution in band. 




f<Hz) 



FIG. 2: Several precession cycles contribute comparably to the 
signal: Power {\hx\^ /Sh{f)) versus frequency for a spinning bi- 
nary viewed edge on in the x polarization, computed directly from a 
Fourier transform of the gravitational-wave signal. In this figure, the 
binary is a 10 + IAMq BH-NS binary with xi = 1 and /3* « n/3 
(solid). Each pair of peaks corresponds to a single precession cycle. 
The signal is divided by the high-power zero-detune advanced LIGO 
noise PSD (power spectral density) B4I1 . 



expression, combined with the leading-order Newtonian ex- 
pression r(/) = M(M7r/)^2/3. 



V 



2.5xi{mi/10m2) 
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1/3 



(4) 



Figure [T] shows an example of /?(/) deduced from this expres- 
sion. 

During simple precession over longer timescales, the width 
of the precession cone increases on the gravitational-radiation 
timescale (rg^u ~ [r''/(77Af'^)](5/64)). During this slow in- 
crease, the orbital angular momentum traces out an ever- 
widening spiral at the precession frequency: 
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where we have replaced r by the leading-order PN expression 



= (M/r) = (Af7r/)-2/3 and Ln = piM{MTrf)-^ 



/3. 

ACST provide an explicit, algebraic solution for the spins as 
a function of time [their Eqs. (59-63)]. 




FIG. 3: Precession opening angle versus spin-orbit misaligne- 
ment: Relation between k — L.S, and fi = cos^^ L.J for 7 = 
ISI/IL] = 0.5, 0.6 ... 1.5. Binaries that are angular-momentum- 
dominated (7 < 1) can have only a small range of precession cone 
opening angles, bounded by Pmax ~ cos"'^ ^1 — 7^. Binaries that 
are spin-dominated (7 > 1) can have all possible precession cone 
opening angles. For extremely spin-dominated binaries n = cos/3. 



Regions of parameter space I: L,S, dominated or 
intermediate 



Single-spin binaries can be loosely subdivided into rare 
transitional precession and ubiquitous simple precession 
[ACST]. In band, simple precessing binaries can be either L 
or S dominated. A spin-dominated binary can have an arbi- 
trary opening angle. An L dominated binary, by contrast, has 
a precession cone opening angle necessarily smaller than tt/2. 




FIG. 4: Angular momentum dominated versus spin-dominated 
binaries: In terms of the mass and spin of the black hole, contours 
of the ratio 7 — \S\/\L\ evaluated at 40 Hz. The bottom left re- 
gion is angular-momentum dominated (\L\ ^ |S|); the top right 
region of large black hole mass and spin is strongly spin-dominated 
i\S\ > \L\). Contours show the ratio |S|/|L| ee 7 = 1 (thick 
curve), sin7r/4 = l/\/2, sin7r/6 = 1/2 and sinvr/S, evaluated 
with a 1.4Mq NS companion, versus the black hole mass and spin 
parameters rriBH, Xbh- Above (below) the thick curve, BH-NS bi- 
naries' total angular momenta are spin (orbit) dominated in band. If 
spin and orbital angular momenta are nearly antialigned, these bi- 
naries have undergone transitional precession at lower frequencies, 
typically not in band. Conversely, for orbital-angular-momentum- 
dominated binaries (7 < 1), transitional precession has not occurred 
in the past at lower frequencies and may, if anti-aligned and 7 near 
1, occur in band in the immediate future. Finally, below the bottom 
curve, BH-NS binaries waveforms are modulated little by precession 
in band. 



and bounded above by 



Pmax = sin ^7 



(8) 



In the neighborhood of this extreme misalignment, at k = —7, 
the opening angle is nearly stationary with spin-orbit mis- 
alignment (i.e., d cos f3/dK ~ 0). In short, a distribution of 
L dominated binaries has two choices for spin-orbit misalig- 
ment (i.e., two values of k) consistent with each realized open- 
ing angle. Additionally, because of the local maximum in /3 
as a function of k, a randomly oriented distribution of spins 
will have opening angles /3 that cluster near that maximum 
(i.e., /3 ~ Pmax)- To illustrate which regions are L and S 
dominated. Figure HI shows contours of constant 7, assuming 

TO2 = 1.4Af0. 



Regions of parameter space II: Steady precession and 
geometry 



Unless transitional precession happens in band, ground- 
based gravitational-wave detectors are sensitive to emission 
from a relatively well-defined epoch: the precession cone has 
relatively constant opening angle [Fig. [Tl. Quantitatively, 
we define a reference frequency fpeak corresponding to the 
frequency up to which half of the signal power has been ac- 
cumulated. The specific reference frequency depends on the 
noise curve adopted.^ For this paper, we adopt the fiducial 
advanced LIGO noise curve with zero-detuned signal recy- 
cling; see 1411). This includes a low-power mode for which 
fpeak ^ 40 Hz and high-power for which fp^ak ^ 60 Hz. 
However, all planned noise curves we have examined have a 
reference frequency in the neighborhood of which a constant 
precession cone is a good approximation. Henceforth the ra- 
tio 7 — |S'i|/|L| and opening angle /? between L and J will 
refer to quantities predicted at th is frequency by the simple 
precession expressions [Eqs. 4 3 1.^ 



Second, not only is the precession cone nearly fixed, but 
as shown in Figure l2] at least a few complete precession cy- 
cles occur between 20 — 100 Hz, where most of the signal-to- 
noise accumulates. For example, for an angular-momentum- 
dominated binary (7 <C 1), the number of precession cycles 
for a single-spin binary can be approximated by the spin- 
independent expression 



Np ~ / dforb-Tj. ^p 

Jirfmin "Jorb 

27(l + 0.75m2/mi) 

M/1OM0 



{M^fmax)-^] 



(9) 



with a comparable but spin-dependent number for an S- 
dominated binary (7 ^ 1); see ACST Eqs. (45, 63) for 
a general solution. As indicated by Figure |2] each preces- 
sion cycle usually accumulates a comparable proportion of 
detectable power (i.e., each pair of peaks is a similar order 
of magnitude in area). More critically, the figure indicates 
that at least one and often several precession cycles contribute 
to the total signal to noise. With many precession cycles, a 
gravitational-wave detector should be relatively insensitive to 
the initial value of the precession phase. 

For our purposes, then, the binary undergoes nearly steady 
simple precession in band. The instantaneous beam pattern of 



^ In the text we choose the reference frequency as the half-power point, 
where J f~'^'^/Sfi{f)df reaches half of its total value. Alterna- 
tively, the reference frequency can be set by maximizing dp/d In / = 
4/|/i(/)p/Sh, or even phenomenologically, in whatever maimer is 
needed for numerically-calculated amplitude and match to reproduce our 
expressions. For the noise curves considered in this paper, all approaches 
nearly agree. 

* For simplicity, we adopt the leading-order (Newtonian) expression for 
r(/). Higher order coiTections are small. 




FIG. 5: Coordinates for steady precession: View from the side 
(top panel) and along the line of sight (bottom panel) of the ge- 
ometry of a typical precessing BH-NS orbit, on timescales t long 
compared with precession but short compared to gravitational-wave 
decay (Qpt 2> 1 3> t/Tgw). The light shaded region indicates the 
cone swept out around the axis by the orbital angular momentum (the 
direction of strong, circularly-polarized gravitational-wave emission 
from each orbit). In both panels, the blue line indicates a direction 
perpendicular to L. The dark shaded region indicates the correspond- 
ing region swept out by the orbital plane (the direction of weakest, 
linearly-polarized emission); observers along these lines of sight see 
strong amplitude modulation. The angle /3 between the total angular 
momentum J and L is nearly constant on precession timescales. On 
longer timescales, as /? increases, the line of sight to the observer 
is often enclosed in one of the two shaded regions. On the bottom 
panel, the x and y axes correspond to the arms of the detector, pro- 
jected into the plane of the sky. The dotted lines correspond to "null 
lines", the directions that, when L is parallel to them, produce zero 
amplitude in our detector; these lines lie in the plane of the sky, per- 
pendicular to our line of sight. [To distinguish J and L in the plane 
of the sky, the two panels adopt different reference times.] 



the binary is aligned with the instantaneous L, and is the same 
beam pattern as a non-spinning binary |42]. Along the axis 
aligned with iL, the radiation is circularly polarized and in 
the orbital plane, the radiation is linearly polarized. Over a 
longer timescale, L precesses around J, sweeping the beam- 
pattern around the precession cone. Any given line of sight 
therefore can be characterized with its proximity to the orbital 
plane (an amplitude minimum and specific, linear polariza- 
tion) and ±i (an amplitude maximum and specific circular 
polarizations). 



D. Constant precession cone coordinates 

In steady precession the unit vector L rotates regularly 
about J. In practice, however, we interpret the gravitational- 
wave strain relative to a frame associated with our line of sight 
— n, where fi is the emission direction from the binary. To sim- 
plify calculations, we adopt a frame defined by this direction: 
z = n, along with two perpendicular directions .x, y. Adopt- 
ing an initial phase so the projection of L into the plane of 
the sky (i.e., perpendicular to n) is along J, the unit vectors 
needed to describe the steady precession approximation are 



J 

P 



sin 9 J [cos i^jx + sin ^jy\ + cos Ojfi 
J X {J X n) 



(10a) 



I J X (J X n)\ 
= cosOjlcostpjx + sinipjy] — sinOjh (10b) 

Here tpj is the orientation of J projected into the plane of the 
sky and 9j is the angle between J and h; see Figure l5] In 
terms of these vectors, the precession of the orbital angular 
momentum about J can be described as 



L 



cos f3J + sin (3R j {a)p 



(10c) 



where Rj{a) is a rotation operator about J and a — J flpdi 
is the precession phase of L around that axis ||35]| . 

The gravitational-wave strain recovered by a detector de- 
pends on the orientation of L relative to the radiation frame. 
In terms of the coordinates above, the inner products needed 
are 



L ■ n 
L ■ X + iL ■ y 



~ cos /3 cos 0j — sin /3 sin 6'j cos a (lOd) 



— J^j 



e^^'' [sin /3(cos a cos Oj + i sin a) 



cos /3 sin 0j 

J4' 



1 - (L • n)2 



(lOe) 
(lOf) 



where ipL is shorthand for the instantaneous orientation of L 
on the plane of the sky. 

For convenience and without loss of generality, in what fol- 
lows we will assume the binary is directly overhead a single 
interferometer, with arms along x and y. For this orientation, 
the detector response functions are F^ = cos2ipj,Fx = 



sin 2^ J. If the source is not directly overhead, the single- 
detector response functions 

F-|_ = - (1 + cos^ 9) cos 2(f) cos 2^p — cos 6 sin 2(f> sin 2-!/) 

(11a) 

Fx — cos 9 sin 20 cos 2tp -\ — ( 1 + cos^ 9) cos 2(f> sin 2')/' 

(lib) 

can be rewritten as an overall scaling and a polarization shift 



Fx =Fosin2(^ + Vo) 



(lie) 
(lid) 



with 



Fo = \/ ((1 + cos2 6l)/2)2 cos2 2(f) + cos2 6* sin^ 20 (1 le) 

cos6' 



tan 2-00 



(l + cos2 6l)/2 



tan 20 



(llf) 



III. CRITICAL VIEWING ORIENTATIONS AND 
DOMAINS 

The instantaneous emission pattern of gravitational radia- 
tion is aligned with the instantaneous orbital angular momen- 
tum L and is the same as the beam pattern of a non-spinning 
binary i.e. an I ~ |to| — 2 quadrupolar configuration. The 
radiation is circularly polarized along L and — L, and lin- 
early polarized in the plane perpendicular to L. In particular, 
when the orbital angular momentum is instantaneously per- 
pendicular to the line of sight, one (linear) polarization of the 
gravitational-wave signal is instantaneously zero. The con- 
dition that one polarization be zero at some instant defines 
a two-dimensional surface in the three-dimensional space 
(9 J, if) J, f3) of all possible precession cone geometries. As we 
show below, this surface decomposes the three-dimensional 
space into three distinct regions, corresponding to different 
ways the orbital plane crosses the line of sight or, equivalently, 
different ways the orbital angular momentum wraps around 
two null lines. 

The detector is not sensitive to (leading-order) gravitational 
radiation emitted from the binary when the orbital angular 
momentum simultaneously satisfies the following two condi- 
tions: 



= L-n^O 

= (L-x)^ -{L- yf 



(12a) 
(12b) 



Solutions to this expressions exist if and only if the angle be- 
tween J and one of the four critical directions indicated in 
Figure |6| is equal to the angle /3 between ,/ and L, or equiva- 
lently if the following expression is zero: 

dnull = (J-Wl -COS^)(J-'yi +C0S/?) 

X (J • i;2 — cos/3)(J • t;2 + cos/?) (13) 

= [{J ■ vxf - cos^ (i][{J ■ V2f - cos" p] (14) 



where vi — [x + y)/^/2 and V2 — {x — y)/^/2. Substi- 
tuting in the coordinate form for J given above and solving 
dnull = for cos2 20; leads to an expression for the boundary 
polarization angle 'if'J,crit between different regions: 



C0S^(2V'J,crJt) 



4cos2^(cos2/3-l) 



sin^e 



(1 + cos 2/3) (cos 2/3 + cos 29) 



(15) 



These configurations defined above divide the space of ori- 
entations {9j,f3,if;j) into five regions, henceforth denoted 
Rl±, R2±, and R3 [Figurel7|. As our notation suggests, us- 
ing symmetry we can relate two pairs of regions by reflection 
symmetry /3 — > /? — tt. The region Rl (the union of Rl-|_ and 
Rl_) corresponds to all configurations such that the orbital 
angular momentum does not wind around either null line. It 
includes all points where the orbital plane does not cross the 
line of sight, plus those connected orientations which never 
have a null in the waveform amplitude. Due to a coordinate 
degeneracy in FigurelT] the region shown corresponding to Rl 
consists of two disjoint mirror-image copies Rl±, related by 
/3 — > TT — /3. These two regions include the special cases of 
spins aligned or antialigned with the total angular momentum, 
viewed along or opposite to the the total angular momentum. 
Conversely, the region R3 corresponds to binaries and view- 
ing geometries so L winds around both null lines. The region 
R3 includes the special case of a precession disk (/3 — tt/2), 
viewed favorably edge on (9j = 7r/2, 0j = 0), plus those 
connected orientations free from amplitude zeros. The two 
regions R2-|- correspond to the ways L can wind around one 
or the other null line. At every polarization except tfjj = 0, 
the regions Rl and R3 are separated by R2. For ipj = 0, the 
three regions overlap at the zeros of 



cos 29 J + 2 cos 2/? + 1 = 2[sin^ 9 + cos 2/3] 



(16) 



IV. WAVEFORMS IN CONSTANT PRECESSION CONE 
LIMIT 

In this section we review how precession modulates the 
amplitude and phase of the binary's gravitational-wave sig- 
nal. We express precession-induced amplitude and phase 
modulations in terms of a complex factor z. We show the 
gravitational-wave phase has both oscillating and secularly 
increasing contributions from precession. The secular effect 
arises as a fraction of the precession phase is imprinted on 
the signal. We demonstrate that each of the previously iden- 
tified regions corresponds to a different amount of secular 
precession contribution. To leading order, gravitational-wave 
detection with non-precessing templates works by match- 
ing a signal with similar phase evolution. In this sense, 
gravitational-wave signals from each of the three regions are 
phenomenologically similar (and can be phenomenologically 
distinguished from one another), as each region has a charac- 
teristic secular phase evolution change from precession. 

Explicit expressions for the leading-order gravitational 
waveforms produced by instantaneously stationary orbits are 




FIG. 6: Three regions have three types of precession cones: The figure above shows the precession cone for three fiducial binaries, char- 
acterized by the angles {6j, j3, tpj) = (7r/4, n/6, 0), (7r/2, 7r/4, 7r/4) and (7r/2, tt/S, 0). Each figure shows the precession-induced path of 
L{t) (colored circles), the total angular momentum direction (arrow), and an example of an instantaneous L direction (arrow). As these three 
examples demonstrate, each region corresponds to a different way the precession cone wraps around the two dotted lines. These two lines 
correspond to orientations of L such that the detector would instantaneously see exactly zero amplitude. We will henceforth denote them as 
"null lines". These three precession cones are representative examples of the three regions discussed in the text: Rl (left panel), R2 (center 
panel) and R3 (right panel) 
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FIG. 7: Region divisions in coordinate form: The three ways the 
binary's orbital angular momentum can wind around the two null 
lines are separated by cases where the cone is tangent to one or the 
other axis. In this figure we express this dividing surface using ex- 
plicit coordinates for the 3d configuration space: z = cos^ 2ipj, 
X — cos26j, y — cos 2/3. Points on this surface are simultaneous 
solutions to Eq. |T2), as calculated by Eq. (Tsl. In this figure, the 
Rl region is the largest, farthest from the viewer; the region R3 is the 
region above the surface, closest to the viewer; and the region R2 is 
below the surface. 



FIG. 8: Phase evolution versus time: For three selected viewing 
geometries of a single binary with /3 — 60°, a plot of $ — $s(i) = 
arg(2) versus a, the precession phase. This plot shows how much of 
the precession phase is accumulates in the gravitational-wave signal 
due to precession [Eq. l|21[(]. These three phase evolutions are gener- 
ated from Eqs. ifTTJ. The dotted line corresponds to Wa, where W 
corresponds to the rate of secular phase increase appropriate to that 
line of sight. The three viewing geometries lie in the three regions 
R1,R2,R3 described in the text. 



available in the literature; see, e.g., Eqs. (31-37) and Fig. 1 
from Buonanno et al. 1201 . Additionally, though these expres- 
sions can be evaluated in any coordinate system (see, e.g.. 
Fig. 1 in Buonanno et al. |20|), we adopt a specific coordi- 
nate system aligned with the total angular momentum. (These 



coordinates are different from the frame adopted in Eq. ( 10 1.) 
In our coordinate system, the expressions provided in Buo- 



nanno et al. |i20J correspond to the expressions below, with 
their l corresponding to our /3 and their 8 corresponding to 
our 6 J. In these expressions, the angle a is the accumulated 
precession phase j flpdt of the orbital angular momentum 
around J. The reference phase <I>5 used in BCV is given by 
4>orb — Oi COS /3 for cf)orb the orbital phase: 



h ex (C+i^+ + Cx-Fx) cos(2(/)orfc — 2a cos /?) + (S'+F-f + S'x-Fx)sin(2(/)or6 — 2a cos/ 
ex Re[(C+Fx +CxFx) -z(5+F+ + 5xFx)] x e2.(0„,„-aco.;3) 



Cx 



(17a) 
(17b) 



C+ = - cos^ 0j[sin^ a — cos^ /?cos^ a] + -(cos^ (3ain^ a — cos^ o;) — - sin^ 6'j sin^ /^ ~ 7 sin 26',/ sin 2/3 cos a (17c) 
S+ = -(1 + cos^ 0j) cos/3sin2Q! + - sin20jsin/3sinQ; 



- cos 9{1 + cos^ /?) sin 2a — - sin 9j sin 2/3 sin a 
' cos 6'j cos /3 cos 2a — sin 6j sin /3 cos a 



(17d) 

(17e) 
(17f) 



r 



with a corresponding expression for the orthogonal polariza- 
tion. For the calculations needed here, however, the waveform 
is most usefully expressed as amplitude and phase modulation 
of a sinusoid: 



h ex Re[Ae2**"'"-=] = Re[ze^''^'f''" 



b — QCOS;3)l 



<1> 



^orb 



acos/3 + -argz 



(18) 
(19) 



z = iC+F+ + C^F^)-i{S+F+ + S^F^) (20) 



Practically speaking, the phase increases and oscillates with 
time. The gravitational-wave phase ^wave can be decom- 
posed into three parts: 



^wave 



^orb 



Wa 



Kprec 



(21) 



These three parts are (a) orbital modulation (porb', (b) a 
precession-induced secular increase in phase imposed by geo- 
metric effects associated with precession (Wa); and (c) resid- 
ual modulations in phase caused by precession {5(t)prec)- Both 
the secular evolution factor W and modulations depend on the 
line of sight; see FigurejHlfor three examples. 

The secular phase evolution W is particularly simple. For 
example, for binaries and viewing orientations where the or- 
bital plane never crosses the line of sight, corresponding to a 
single helicity being present along the line of sight, the secular 
phase must accumulate steadily in proportion to this helicity. 
Using the special case 6* = 0, tt to determine suitable coeffi- 
cients, one can show that 



W 



(sign(cos/3) — cos/3) 



(22) 



More generally, the secular phase factor 2W is —2 cos /3 plus 
an integer, related to the number of times the complex number 
z winds around the point z = in a precession cycle (i.e., 
versus a): 



W = -n^tnd 



cos P — — cos j3 



1 

A'Kl 



dz 



(23) 



where this integral follows a contour in z corresponding to a 
evolving from to 27r (i.e., one precession cycle). 

Figure [9] show examples of different trajectories for the 
complex amplitude z{t). Depending on how the angular mo- 
mentum winds around the two "null lines" [Figure l6), the 



complex number z derived from the angular momentum and 
the detector orientation can wind around the origin (the point 
corresponding to L along those null lines). This integral can 
be explicitly evaluated, for example, by substituting u — e'" 
and performing a contour integral. Since by construction the 
complex number z is exactly zero at some point on its cycle if 
and only if Eq. ( 12 1 holds, the connected subsets of Rl, R2, 
and R3 are also regions of constant winding number n^wind- 
Keeping in mind the regions have different connected subsets 
depending on the sign of cos /3, we evaluate W and find 



W 



■ cos /3 



sign(cos/3) 
^ sign (cos /3) 




Rl 
R2 
R3 



(24) 



ACST find an approximately comparable result [their Eq. 
(65); see also and |37|] by integrating an ODE for the phase 
perturbation due to precession. That expression does not in- 
clude the region R2 or correctly identify the conditions that 
cause transition from one secular phase to another. 

The complex amplitude z{t) also encodes the way that the 
gravitational-wave phase and amplitude "oscillate" about their 
representative values. Binaries whose angular momenta are 
aligned with their orbit do not precess, so they do not accumu- 
late secular precession phase (W = 0, as z ex e*^") and their 
phase does not oscillate. As is apparent from the z{t) trajec- 
tories in Figure [9J any oscillations in amplitude and phase are 
more extreme for trajectories passing closer to z ~ 0. These 
extreme oscillations cannot be well-fit by a non-precessing 
template. We will show below that the best fitting factor be- 
tween a precessing signal and a non-precessing template fam- 
ily correlates directly with the amount of phase and amplitude 
oscillation a signal exhibits. 

By way of example. Figure l8] shows how precession intro- 
duces additional secular and oscillating contributions to the 
gravitational wave phase of the same binary, seen along dif- 
ferent lines of sight. The three curves in this figure show 
arg(z) versus a, as extracted from z(a) [Figure l9). We can 
clearly see the distinct secular phase trajectories, correspond- 
ing to each of the three ways L can wind around the "null 
lines" (equivalently, how z{t) can wind around the origin). 
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FIG. 9: Complex amplitude contours: Three examples of the complex amplitude z{a) over a precession cycle, in the constant precession 
cone approximation. The three configurations shown correspond to the three fiducial precession cones shown in Figure[6] 



V. AVERAGING AMPLITUDE AND MISMATCH WITH 
THE PRECESSION CONE 

In this section we estimate the response of gravitational- 
wave detectors and non-precessing data analysis strategies to 
a precessing signal. Using a separation-of-timescales argu- 
ment, we show that precession-induced modulations decou- 
ple from the orbital phase, allowing key expressions like the 
signal amplitude and mismatch against proposed signal tem- 
plates to be computed by an average over the precession cone. 
We apply this idea to compute how precession modifies two 
quantities: the signal strength (SNR) and the mismatch of 
a non-precessing template with a precessing signal. Com- 
bined, these expressions tell us the relative signal strength a 
data analysis pipeline would recover from a precessing binary, 
given we know how well it performs for a non-precessing 
counterpart (i.e., one with similar masses but aligned spins). 

The signal-to-noise ratio for the detection of an exact tem- 
plate waveform h^ using an optimal filter constructed from a 
model waveform hm used as a template is given by 



Pm = {he\hm) = 



{he\hm) 



(25) 



where the noise-weighted inner product {he\hm) is given by 



{h,\h^) = 2 / df 



Mf)h*mif) 
Sn{\f\) ■ 



(26) 



Here i^nd/j) is the one-sided power spectral density of the 
detector strain noise n(t), defined by 

{n{f)n*{n) = ^6{f-f)S„{\f\), (27) 

and h{f) is the frequency-domain representation of the wave- 
form, given by 



Hf)^ / /i(i)e-'"^* di 



(28) 



If we know the exact waveform, the strength of the signal 
is given by p^ — {he\he). However, if there is an error in the 



waveform model h„ 
by a factor 



, then the signal-to-noise ratio is reduced 



p,n = Mp 



(29) 



where < 7\f < 1 is the match of the model signal. The 
model waveform /i,,, — /i„i(A) is a function of the model's 
intrinsic parameters A (e.g. masses, spins, etc.) and extrinsic 
parameters Ag (e.g., time of arrival, coalescence phase, sky 
location, etc.). The model waveform having the same physi- 
cal parameters as the exact waveform ft,,,, might not have the 
highest match, and so we define the fitting factor F as 



F 



niax(fte 

A.Ae 



^m(A, Ae)). 



(30) 



For the best-fit parameters, the signal-to-noise ratio will be 
reduced by a factor of pm = Fp. 

In this section we use the constant precession cone ap- 
proximation to derive closed-form expressions for the fitting 
factor F and the geometrical factors in the signal amplitude 
p^. Formally, our approximation arises through separation of 
timescales. At each instant, the gravitational-wave emission is 
instantaneously quadrupolar Expanding in powers of ilp/irf, 
the time-domain amplitude modulation A{t) implied by the 
beam pattern translates directly to frequency-domain modula- 
tion A{t{f)) [ACST Eq. (36-37)]; see, e.g.. Figure [I] Each 
frequency identifies an orientation of the angular momentum 
vector L{f) and thus precession phase a. So long as a few 
precession cycles dominate peak emission, a single (effec- 
tively "constant") precession cone dominates as well. Adopt- 
ing this precession cone as fiducial, we can transform the inner 
products involving A{t{f)) — Ao{f)B{L{t)) into inner prod- 
ucts involving the "slow" degrees of freedom (i.e., just A{f)) 
and averages over the "fast" (precession) degrees of freedom 
{B{L)). The averages over slow degrees of freedom do not 
explicitly depend on precession geometry, detector noise, or 
component masses. The expression factors. We therefore ar- 
rive at expressions for p'^ and F in terms of averages. As these 
averages encode all orientation-dependent effects, we empha- 
size the averages below, keeping in mind that their prefactors 
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FIG. 10: Amplitude ratio model: For spinning precessing binaries, the function s{9j, I3,^j — 0) describing the average amplitude along 
a given line of sight, assuming a steadily precessing binary with precession cone opening angle /3. The function s{6j, I3,ipj) [Eq. pTb] 
is approximately the ratio p/ pns between the amplitude p of a spinning, precessing binary seen in-band in a configuration and line of sight 
specified by {6j, P, tpj) and the amplitude pns of a nonspinning binary of comparable masses. 



(if needed) can be easily computed once and for all, using a face-on binary. 



A. Averaged amplitude 

The precession of the beam across the line of sight modulates the amplitude seen at a single detector, in both sin and cos 
quadratures. We estimate the signal power by adding each quadrature's contribution independently, including the prefactor 
implied by the orientation of L relative to the line of sight and detector. We then average the resulting expression over the 
entire precession cone. Compared to the geometrical factor for an optimally oriented non-precessing source, the amplitude p^ is 
smaller by a factor s^: 

Ae„p,^,) . / (^ + ^^^ ■ -^^)^ cos- 2Mt) + (L ■ nr sin- 2^,(t)\ J ^^ + ^^Z ^^'^' - ^' ' ^^^ ' ^^'^' sin^ 2V.. (t)^ 



{l + iL-nf) 



2\2 



{L-xf{L-yf 



1024 



[{cp{x - if + x-}(35y- + lOy - 13) + 2x (5t/ + 166?/ + 53) - 13y- + 106?/ + 45l] 



where a: — cos26j,y = cos2/3 and Cp — cos-iipj. Note that all polarization dependence enters through a term oc Cp(l— a;)-+x- 
the average amplitude is almost polarization-independent near dj ~ 0. 



Figure 10 shows how the relative amplitude of a BH-NS binary seen directly overhead a single detector should change, versus 
the purely geometrical orientation parameters 9, /3,ipj describing the constant precession cone. 



B. Averaging and fitting factor 



The phase evolution is more subtle. As simple preces- 
sion parallel-transports the orbit along a precessing path, the 
gravitational-wave phase (j)wave of the single polarization to 
which our detector is sensitive oscillates about a combination 
of the orbital phase and a secularly accumulating proportion 



of the precession phase a: 

'24'wave = '^4>orb + '^Hprecio) + 2W a 



(33) 



as described above. Finally, the spin also influences the accu- 
mulation of orbital phase (t>orb'- the projection of Si along L 
produces a "gravitomagnetic repulsion", changing the orbital 
evolution (porb by slowing the rate of increase in frequency; 
see, e.g., Eq. (6) in |,19J . 
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For detectors like LIGO, which are sensitive to an epoch 
where the binary's opening angle is nearly constant, the sec- 
ular phase changes {(fiorb and Wa) can be well-fit by a suffi- 
ciently generic non-precessing waveform.^ For example, the 
secular phase evolution of many precessing waveforms can 
be well-fit by a non-precessing waveform with different chirp 
mass and mass ratio (M.ci'n)- On the other hand, amplitude 
and phase oscillations in A{t) and 5(j)prec cannot be fit with 
non-precessing templates. 

Using the ansatz that the best non-precessing template ft,,„ 
will perfectly reproduce the secular phase, the time-averaged 
overlap between the spinning waveform he and the best non- 
spinning template must be proportional to the best nonspin- 
ning fitting factor. We therefore define approximation to the 
fitting factor F that captures the effects of precession by 



sT: 



{K{t)\h^[t)) 



(34) 



where s is the average amplitude from Eq. pT|) and where 
the average is over a single precession cycle of L in he, and 
where hm is the best-fitting nonspinning template to the entire 
cycle. Replacing a time average by an average over the pre- 
cession phase a and the true signal h^ (t) with the nonspinning 
template times the amplitude and (fourier domain) precession 
phase terms A and S<j> described above, this expression re- 
duces to an average of the oscillating amplitude and phase 
over the cycle: 



sJ- — max 



00 



^cos(25^ + (/)o)da/(27r) (35) 



We will provide a closed-form expression in the next sec- 
tion. However, to understand that result, it is helpful to think 
directly in terms of this integral. For example, in the limit 
of small amplitude and phase modulations occurring periodi- 
cally on the precession cone, a non-precessing template model 
should fit them to no better than 



J" 



1- 



1 



{{2Scf,r) + {{SA/Ar)-{dA/Ay 



(36) 



[This expression can be rederived directly from the original 
integral.] In particular, the amplitude and phase modulations 
introduced by precession inevitably diminish the ability of a 
non-precessing template to match them, by an amount pro- 
portional to the standard deviation of oscillations in time and 
phase about the reference model. 

Oscillations in phase and amplitude are most extreme for 
binary geometries near the boundaries between regions. A 
non-precessing signal should be maximally unable to match 
the large phase and amplitude variations near this surface. For 
example, a signal on the surface has different epochs, sepa- 
rated by zeros of the (complex) amplitude. A non-precessing 



signal can coherently reproduce the phase trend in one portion 
of the signal or the other, but not both. 



C. Constant Precession Cone 

In the limit of a constant precession cone, the integral 
can be performed analytically. Replacing a time average by 
an average over the precession phase a and the true signal 
h{L) with the nonspinning template times the amplitude and 
(fourier domain) precession phase terms A and Sep described 
above, the angle-averaged expression reduces 



sJ- = y'isTs)^ + {sTeY 



(37) 



where in the second line we convert maximization of the inte- 
grand over (/)o into two integrals. 



sTc 



Atot cos 2S(j) 
Atot sin 264) 



da 

2tt 
da 

2^ 



(38) 
(39) 



This expression reduces to a simple contour integral: 

, (z/Z + {z/Z)*)da 



2 27r 

{{z/Z) - {z/Z)*)da 
2i 2tt 

z{a) ^ A„,e'2(*-*+"^) 
Z{a) 



_ i2Wa 



(40) 

(41) 

(42) 
(43) 



where z is the complex representation of the waveform ampli- 



tude and phase given in Eq. ( 20 1, while dividing by Z is equiv- 
alent to subtracting the secular phase. Using the complex vari- 
able u = e'" to represent Z — v?^ — u"™"><i and the sinu- 
soids in C+.x and S'+^x analytically via cosa = (u+ \/u)/2 
and similarly, we find the fitting factor predicted can be calcu- 
lated using the contour integral 



sF = 
I = 



zjZda 
27r 



z/m" 



^du 



2Triu 



(44) 
(45) 



This trivial contour integral corresponds to identifying terms 
ex m" in z for different n. The value of this contour integral 
depends discontinuously on the winding number riwind- Us- 
ing the definitions of C, S to perform this integral, we find 
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T 



- 1 cos 2^ J sin^ /3 sin^ 9 

(2 sin /3±sin 2/3) (cos 2ipj sin 2e=F2i sin 9 sin 21/J.7) 



= 

= ±1 



(46a) 



{l±cosl3y 



g [cos2'0j(l + cos^ 6*) =F 2icos6'sin2?/'j] n^ind = ^2 

Finally, we estimate the fitting factor for all angles by dividing this closed-form expression by the closed-form expression for s: 

F{e.j,p,^j) = |/(0j,/3>j)|/s(0j,/?,Vj) (46b) 

I 



This ratio varies strongly in all three parameters and is best 
summarized by the explicit expression above. Though com- 
plicated, some important physical lessons can be drawn. First 
and foremost, for non-precessing signals the predicted fitting 
factor is F = 1. Second and most critically, while the pre- 
dicted fitting factor has a local extremum with J^ ~ 1 in each 
of the three regions, the predicted fitting factor drops rapidly 
as we approach the surface dividing them: at the surface de- 
fined by Eq. (15 1 and shown in Figure [7] In other words, 
this expression suggests the worst possible matches will oc- 
cur when the angular momentum at some point in the preces- 
sion cycle is along one of the two null lines. As noted above, 
the low fitting factor in this region occurs because the signal 
consists of comparable amounts of power with two competing 
and incompatible secular phase trends. Finally, under special 
conditions such as 6' ~ in R2 and R3, this expression pre- 
dicts pathologically low matches P ~ 0. As we will discuss 
at greater length in a forthcoming publication, an empirically 
more successful approximation replaces the options in / by 
the the maximum over the three options. 

Though derived under the seemingly-strong expression of 
a constant precession cone, we have verified this expression 
applies to a broad family of real BH-NS inspirals, by calcu- 
lating the overlap between the gravitational radiation emitted 
from non-precessing and precessing binaries and maximizing 
over the non-precessing signal's parameters. We will describe 
these Monte Carlo studies in more detail in a subsequent pub- 
lication. In brief, these studies corroborate many earlier stud- 
ies that examined how well non-precessing signals recover 
precessing waveforms ll20l ISTl l37l . First, in special but not 
uncommon circumstances, fitting factors as low as 0.6 occur. 
Second, small fitting factors (i.e., large mismatches) can occur 
even with small amounts of spin-orbit misalignment, though 
they become more common as misalignment increases. For 
nearly aligned binaries, signals with the lowest match corre- 
spond to directions with low signal power. Contrary to pre- 
vious claims, however, we identify some configurations with 
both low fitting factor and significant signal power along the 
line of sight. 

This analytic study differs from all previous investigations 
in that we assume a non-precessing signal model can always 
match any secular phase evolution; identify which precessing 
configurations lead to large modulations; and quantify the ex- 
pected purely geometrical effect these modulations have on 




FIG. 11: Distribution of fitting factors for single-spin BH-NS: 

Using only the analytic expressions provided in this paper, the pre- 
dicted fraction P{> J-) of fitting factors greater than J^ for a ran- 
domly oriented BH-NS binary with mi — lOM0,m2 = IAMq, 
with the dimensionless black hole spin xi randomly oriented and 
with magnitude uniformly distributed between and 0.98. Compare 
to the solid line in the bottom right panel of Figure 9 in 1311 , although 
there the neutron star also has a spin up to 0.3. 



search is applied to randomly-oriented BH-NS binary. De- 
spite the simplicity of our model, we recover quantitatively 
similar results to previous studies |31 1. Unlike earlier studies, 
however, our expressions allow us both to identify precisely 
which binaries are least well fit and to understand precisely 
why. 



D. Validating our approximations with Monte Carlo 

To demonstrate that our approximation works, we have 
compared randomly-chosen BH-NS binaries with a standard 
non-precessing template bank: Taylor F2 templates with 
physical parameters, hexagonal template bank chosen for 97% 
minimal match; and a likely early advanced LIGO noise 
curve, the low-power zero-detune PSD BTI such that fpeak — 
40 Hz. 



the fitting factor. As an example, in Figure 11 we predict the 
distribution of fitting factors expected when a non-precessing 



For example. Figure 12 shows the results for 30000 BH- 
NS binaries (IOM0 + I.4M0) with random orientations for 
K — L ■ S > —1/2 and spin magnitudes ai £ [0, 1]. The 
limit on k removes some systems with transitional preces- 
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sion for clarity; in the next figure, we allow nearly the full 
range of k. The bottom panel shows the configurations of the 
worst- fit sources for two polarizations (J nearly aligned with 
the detector on the plane of the sky and 45° off the detector 
arms). Both cluster near the predictions of our model (dot- 
ted lines): our approximation correctly identifies the worst-fit 
locations. More generally, the top panel shows our estimate 
works (nearly) everywhere: it compares our fitting factor es- 
timate T to the fitting factor F calculated by comparing that 
signal the template bank. Despite known systematic differ- 
ences between the signal and template model even for aligned 
systems, our match prediction works well almost everywhere. 
For comparison, the vertical dotted line at —0.03 shows the 
typical mismatch expected from template bank discreteness. 

In this comparison, we have artificially eliminated the small 
fraction of worst-performing configurations, corresponding to 
sources recovered with templates near the edge of the tem- 
plate bank where 77 = 1/4. The secular phasing due to the 
aligned component of the spin is best fit by templates with r] 
biased to larger than its true value. Due to the finite parameter 
extent in ry of physical Taylor F2 templates, the conventional 
signal model fails to capture even the secular effects of large 
aligned spins. An aligned but non-precessing template bank 
would however be able to caputre these signals. To faithfully 
reproduce signals with aligned spins using standard template 
banks, we must extend the bank into 77 > 1/4. This nearly 
captures the secular effect of aligned spins, and is a reasonably 
good subsitute for a true aligned spin template bank. Figure 
\T3\ shows the results of this extended-77 bank, compared to a 
much larger set of BH-NS binaries. The top panel shows that 
our approximation continues to work well for generic sources, 
even given systematic differences between signal and tem- 
plate and the absence of spin. In particular, the bottom panel 
demonstrates that except for rare outliers, our approximation 
correctly identifies the worst-fit sources. 



E. Beyond the constant precession cone 

Though we emphasize the value of a (nearly-)constant 
precession cone in approximating binary inspiral, the con- 
cepts and expressions we provide apply equally well to many 
simply-precessing binaries. As an example, in Figure 14 



we 

show the precession-induced phase evolution for two systems 
with slightly different f3, such that they are near a transition 
between regions. One is in Rl and the other in R2. The near- 
ness to the transition surface causes large oscillations around 
the secular increase in phase. Binaries whose opening angles 
/3(/) evolve across these boundaries and particularly across 
tt/2 in band can undergo dramatic transitions in the properties 
of their observable waveform, accumulating only one helicity 
early on and another helicity later. 
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FIG. 12: Comparison of the analytic prediction of Sec. V.C with 
numerical results I: Single binary masses and spin magnitude Top 
pa«e/: For 30000 BH-NS binaries (lOAf© -1-1. 4Af0 with ai G [0,1]) 
with random orientations for L ■ 5 > —1/2, a comparison between 
our estimate J^ and the fitting factor F calculated by comparing that 
signal with a standard single-detector template family (Taylor F2 
templates with physical parameters; hexagonal template bank cho- 
sen for 97% minimal match; early-stage advanced LIGO noise curve 
such that fpeak — 40 Hz 141]). For comparison, the vertical dotted 
line at —0.03 shows the typical mismatch expected from template 
bank discreteness. Despite known systematic differences between 
the signal and template model even for aligned systems, our match 
prediction works well almost everywhere. To highlight its effective- 
ness, we have eliminated the small fraction of worst-performing con- 
figurations, corresponding to sources recovered with templates with 
0.24 < rj < 1/4. Due to the finite parameter extent of physical Tay- 
lor F2 templates, the conventional signal model cannot reproduce 
large aligned spins. Bottom panel: For the same comparison of tem- 
plate bank against signals as above, the configurations of the worst- 
fit sources {F < 0.7) with cos 4?/; > 0.9 (blue; J nearly aligned 
with the detector on the plane of the sky) and cos 4^ < —0.9 (red; J 
nearly 45° off the detector arms, in the plane of the sky). For compar- 
ison, the blue and red solid lines are the surfaces cos 2/3 + cos^ 6 = 
(blue) and cos 29 j + 2 cos 2/3 (red) [Eq. lfT6|]. These lines are 
two one-dimensional cuts through the surface separating Rl and R2, 
where our expression predicts the worst single-detector matches will 
occur. 



VI. CONCLUSIONS 

In this paper we have explored two features of the 
gravitational-wave signal seen by ground-based detectors in 



a "constant precession cone" limit. We compare their ampli- 
tude to similar face-on non-precessing binary; we compare 
their waveform to the best-fitting non-precessing signal. We 
find closed-form expressions for both quantities. Our expres- 
sions involve only geometrical parameters - the orientation of 
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influence that the precession cone has on signal amplitude. 



FIG. 13: Comparison of the analytic prediction of Sec. V.C with 
numerical results II: Monte Carlo. As Figure[T2] except (a) the in- 
jected signal can have generic BH-NS-like masses mi G [3, 15]Mq, 
rn,2 £ [l,5]Af0, dimensionless spin magnitudes x ^ [OjI]' ^^d 
spin-orbit misalignments L-S > —0.9; (b) to better recover aligned- 
spin signals, the template family includes Taylor F2 waveforms with 
77 > 1/4, spaced in chirp mass and -q according to the 2PN hexagonal 
template bank. Top panel: Our match prediction works for generic 
sources. [Because the mass and spin distribution used in this analysis 
favors low 7 and thus binaries nearly aligned in band, this distribu- 
tion partially reflects our ability to recover nearly-aligned-spin bina- 
ries with an extended rj Taylor F2 bank. That said, subsets of the sam- 
ple bounded below in 7 also demonstrate our prediction works well 
for generic binaries with significant in-band misalignment.] Bottom 
panel: Our expression identifies the dominant physical mechanisms 
that produce a poor fit between a spinning binary and a candidate 
non-precessing waveform. Like Figure [T2] the dotted lines show 
the seperatricies between regions. The solid lines show the contours 
predicted from our expression for J^ = 0.75 and -i/i = (blue) or 
-i/i = 7r/4 (red). For comparison, the points indicate all signal versus 
bank fitting factors F < 0.66 for cos 4tp > 0.975 (blue) or F < 0.7 
for cos 4?/) < —0.975. 



the precession cone relative to the line of sight - with no ex- 
plicit reference to the (post-Newtonian) phase model, detector 
noise curve, or binary parameters. In a subsequent paper we 
will use these factors to interpret the selection biases of opti- 
mal and non-precessing searches for precessing binaries. This 
publication will provide explicit Monte Carlo over all possi- 
ble source parameters and sky locations and calculate overlaps 
using state-of-the-art signal models. 

In this paper we provide only the leading-order geometric 




FIG. 14: Discrete phase change: With a fixed viewing geometry, 
the geometrical phase difference arg(2:) on either side of a transition 
surface (here, Rl to R2). As in Figure [8] an entire precession cycle 
is shown. 



Additionally, the same spin misalignments that allow the bi- 
nary's orbit to precess also diminish the binary's inspiral time, 
compared to an aligned non-precessing signal. We will de- 
scribe a more complete model for the spin-dependent signal 
amplitude p^ in a future paper. 

Finally, in this paper we point out that non-precessing 
signals often fit precessing signals' secular phase evolution. 
The secular evolution of (single-detector) gravitational-wave 
phase depends on the line of sight to the source and orien- 
tation of the detector (W). As a result, two independent de- 
tectors running two independent searches may identify two 
very different best-fitting (non-precessing) signals, due only 
to the detectors' orientations. These biases in recovered pa- 
rameters must be taken into careful account when constructing 
coincidence-based search algorithms. We will also address bi- 
ases in recovered parameters in a future paper. 
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